Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 17 February 2008 (MN style file v2.2) 



The real-space clustering of luminous red galaxies around 
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ABSTRACT 

We measure the clustering of a sample of photometrically selected luminous red galax- 
ies around a low redshift (0.2 < z < 0.6) sample of quasars selected from the Sloan 
Digital Sky Survey Data Release 5. We make use of a new statistical estimator to 
obtain precise measurements of the LRG auto-correlations and constrain halo occu- 
pation distributions for them. These are used to generate mock catalogs which aid in 
interpreting our quasar-LRG cross correlation measurements. The cross correlation is 
well described by a power law with slope 1.8 ± 0.1 and tq = 6 ± 0.5/i~^Mpc, con- 
sistent with observed galaxy correlation functions. We find no evidence for 'excess' 
cl ustering on 0. 1 Mpc scale s and de monstr ate that th i s is co nsistent with the results 
of I Serb er et al.l (p006 ) and IStrand. Brunner fc Mverd (|2007[ ). when one accounts for 
several subtleties in the interpretation of their measurements. Combining the quasar- 
LRG cross correlation with the LRG auto-correlations, we determine a large-scale 
quasar bias feggo = 1.09 ±0.15 at a median redshift of 0.43, with no observed redshift 
or luminosity evolution. This corresponds to a mean halo mass (M) ^ 10^^ h~^MQ, 
Eddington ratios from 0.01 to 1 and lifetimes less than 10^ yr. Using simple models of 
halo occupation, these correspond to a number density of quasar hosts greater than 
10~^ft.^Mpc~^ and stellar masses less than 10^^h~^ Mq. The small-scale clustering sig- 
nal can be interpreted with the aid of our mock LRG catalogs, and depends on the 
manner in which quasars inhabit halos. We find that our small scale measurements 
are inconsistent with quasar positions being randomly subsampled from halo centers 
above a mass threshold, requiring a satellite fraction > 25 per cent. 



1 INTRODUCTION 



Quasars are among the most luminous astrophysical objects, 
and are believed to be powered by accretion onto super - 
massive black holes fe.g. lsalpeteilll964l : lL7nden-Belllll969l ). 
They have become a key element in our current paradigm 
of galaxy evolution - essentially all spheroidal systems at 
prese nt harbor massive black holes ([Kormcndv & Richstone 
1 19951 ). the masses of which are correlated with many proper- 
ties of their host systems. However, the physical mechanisms 
that trigger and fuel quasars are still unknown; furthermore, 
it is possible that very different mechanisms dominate at low 
and high redshifts and at high and low luminosities. 

Deep imaging with the Hubble Space Telescope (e.g . 
Bahcall et aLl 'l997': 'M cLure et aLlll999l : lD"unlop et al.ll20oi 



Flovd et air 2004: Letawe et al. 20081 ) suggest that low- 



QSOs are exclusively hosted by bright galaxies with L > L*. 
Radio-loud QSOs reside in early-type galaxies, while their 
radio-quiet counterparts have both early and late type hosts. 



with the fraction of early type hosts increasing with the 
quasar optical luminosity. In most cases there is strong ob- 
servational evidence for th e presence of a youn g (subdom- 
inant) stellar population (jSanchez et al.l [200J) and large 
amounts of gas, irrespectiv e of the morphological type of 
the host (|Letawe et al.|[2007l ). Furthermore, 30-50 per cent of 
quasars appear to be associated with interactions, although 
the number of such imaged systems is small, and signatures 
of mergers are notoriously difficult to observe. The emerging 
picture is that Q SO activity and star formation are inextri - 
cably linked (e.g. [Silverman et al.|[2007l : iNandra et al.ll2007l ) 
in galaxies that contain a massive bulge (and thus a massive 
black hole) and a gas reservoir. 

The clustering of quasars as a function of redshift and 
luminosity provides a different perspective on the above pic- 
ture. The amplitude of clustering on large scales is related to 
the masses of the dark matter halos which host the quasars 
(their environment), which together with the observed num- 
ber density allows us to constrain the quasar lifetimes or 
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duty cycles. The small-scale clustering of quasars can shed 
light on the triggering mechanism for quasars, and the na- 
ture of quasar progenitors. 

However it is only recently that samples of quasars 
have grown big enough (in terms of the number of 
objec ts) to study their clustering with some preci- 



Mvers et al.l l2007al lbl: IShen et al] l2007l : Ida Angela et all 
20081 ). One of the major problems with measuring the clus- 



sion (jPorciani. Magliocchetti fc Norberd 20041: Groom et al.l 
20051: iPorciani fc Norberd l2006l: iHennawi et all l2006l: 



tering of quasars is that they are extremely rare [n ~ 
10"*^ /i^ Mpc"^ at z ~ 0.5). Shot -noise from Poisson fluc- 
tuations in the counts of objects thus obscures their clus- 
tering signal. At low redshifts, this problem is exacerbated, 
requiring measurements in very broad redshift intervals. 



To avoid th is, we cross-correlate a pproximately 2, 500 
low z quasars l|Schneider et al.l |2007|) with a sample of 
450, 00 luminous red galaxies (LRGs) (|Padmanabhan et al.l 
I2OO7I). both selected from the Sloan Digital Sky Survey 
I York et al. I lioOQ . SDSS) and with overlapping redshift dis- 
tributi ons. The LRGs have very reliable photometric red- 
shifts (jPadmanabhan et al.|[2005l ). trace the matter distri- 
bution in a way that is well understood and have a much 
higher volume density (n ~ 10~* Mpc~^) than the quasar 
sample. The cross-correlation can thus be well measured and 
inverted, using the known redshift distribution, to the un- 
derlying 3D clustering. 



In this paper, we make use of several novel 
techniques for measuring the clustering of galaxies 
and quasars, and compute full covariance matrices 
for our estimators from the data themselves. While 
the idea of enhancing the cluster ing signal by us- 
ing cross-correlations is no t new llGroom et al.1 |2004 



Adelberger fc Steidcl 2005a b; Serber et al. 2006; Goil et al 



Strand. Brunner fc Mverj 20071 : Mountrichas et al 



20081 ). the sample size and ability to perform such detailed 
statistical analyses is new to this paper. In addition, the pre- 
cise measurements of LRG clustering allow us for the first 
time to discuss the manner in which both LRGs and quasars 
inhabit dark matter halos at 2: ~ 0.5. 

The outline of the paper is as follows. In !j2]we describe 
the LRG and quasar samples, drawn from the SDSS, that we 
use. The clustering measurements are described in ^ where 
we pay special attention to the techniques used and the er- 
ror estimates. The implications of our results for quasars is 
explored in ^ including comparisons with earlier work. In 
particular we investigate the manner in which quasars in- 
habit dark matter halos at z ~ 0.5. We conclude in ^ Ap- 
pendix|X]contains the technical details of the halo model fits 
used in t his paper, while Ap pen dix [B] recasts the measure- 
ments of ISerber et al.l ()2006l ) and lStrand. Brunner fc Mverj 
l|2007l ) into the framework of this paper, highlighting un- 
appreciated subtleties in their interpretation. Where nec- 
essary we shall assume a AGDM cosmological model with 
fimat = 0.25, = 0.75, erg = 0.8. Also, unless the h depen- 
dence is explicitly specified, we assume h = 0.7. 



2 DATA 
2.1 Quasars 

We use quasars selected from the fourth edition of 
the Sloan Digital Sk y Survey (SDSS) quasar catalog 
(jSchneider et al.l l2007h . This catalog consists of spectro- 
scopically identifie d quasars in t he fifth SDSS data release 
(jAdelman-M cCart hv et al.l [20o3), with an absolute PSF 
magnitude in the i band. Mi < —22.0 and at least one emis- 
sion line with FWHM larger than 1000 kms~^. It does not 
contain Type 2 QSOs, Seyferts or BL Lac obj ects. In order to 
construct a homogeneous sample, we follow iRichards et al.l 
(2006 ) and select objects 

• that were targeted for science 
(SCIENGEPRIMARY = 1). 

• classified by the SDSS photometric pipeline as primary 
(PRIMARY = 1). 

• morphologically consistent with being point sources 
(MORPHOLOGY = 0). 

• with dust and emission line K-corrected i-band magni- 
tudes, i < 19.1. 

In order to cross-correlate with the LRG sample described 
below, we restrict ourselves to quasars that lie within the 
LRG angular mask, and with redshifts between 0.25 < 
2 < 0.6. The resulting sample (denoted ALL below) has 
2476 quasars. One subtlety with the SDSS quasar samples 
are the changes to t he quasar target selection algorithm 
(jRichards et al.ll2002l ) over the lifetime of the survey. How- 
ever, these changes were made to optimize target selection 
at high redshifts, and have little effect on our sample. As we 
discuss below, restricting to quasars selected with the final 
version of the QSO target selection (v3_l_0) does not affect 
any of our results. We therefore do not make a cut based on 
the target selection algorithm. The redshift distribution is 
shown in Fig. [1] while the angular distribution is in Fig. [J] 
Fig. El plots the conditional magnitude distribution of 
our sample. We assume a redshift evolution of M,,i (k- 
corrected to z = 0) given by 



M.,i(2) = M.,i,o - 2.5(fci2 + k2Z 



(1) 



with M^.i.o = -21 .678, fci = 1.39 and = -0.29 
dRichards et al.ll2005l ). where we have converted from the 
SDSS g to i band with g-i = 0.068 ([Richards et al.ll2006l ). 
This defines a sample of quasars (denoted LSTAR) with 
Mi < Mt^i, yielding an approximately volume limited sam- 
ple over the redshift range we consider. We can estimate 
the number density of these samples by integrating the 
broken power-law fit to the QSO luminosity function from 
[Richards et"aLl l|2005l ). 



lQO.i(a + l)(Mi-M,_i) _^ ■]^q0.4(/3+1){A/ — M. i) 



(2) 



If we assume the parameters estimated from the 2QZ and 
2SLAQ surveys i.e. M*,; as defined above, a = —3.31, 
= -1.45, and <! ■* = 5.33 x lO"*' fe^Mpc'Vag*^ 
("Ric hards et^ |2005| ). and no scatter between g and i 
magnitudes, we estimate a number density of 1.7 x 10~^ 
/t'^Mp c"'^. If one adopts the parameters from iBovle et al.l 
(|2000l ). we find a number density of ~ 1.3 x 10"^ /i^Mpc"^, 
approximately 20 per cent lower. Note that we keep the sam- 
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Figure 1. The redshift distribution of tlie quasars [top] and LRGs 
[bottom] used in tliis analysis. The LRG redshift distributions are 
derived from the observed photometric redshift distribution, af- 
ter deconvolving the redshift errors. The dotted lines show the 
redshift distribution for the six individual d^photo = 0.05 LRG 
samples. The vertical lines mark the boundaries of the three 
quasar redshift slices we consider — 0.25 < z < 0.35 (dotted), 
0.33 < 2: < 50 (short-dashed) and 0.45 < z < 0.6 (long-dashed). 



Sample 


Redshift 


^qso 


ALLO 


0.25 < z < 0.35 


435 


ALLl 


0.33 < z < 0.50 


1277 


ALL2 


0.45 < z < 0.60 


1269 


LSTARO 
LSTARl 
LSTAR2 


0.25 < 2 < 0.35 
0.33 < 2 < 0.50 
0.45 < 2 < 0.60 


212 
751 
1094 



Table 1. A summary of the quasar we consider. The columns list 
the sample name, redshift range, and the number of quasars. 



pie definition the same for both these cases, so one is inte- 
grating over magnitudes less than Af, defined by Eq. [1] 

We e stimate the bolomet ric luminosity using the rela- 
tion from ICroom et all l|2005h . 



Figure 2. The angular distribution of our quasar sample, plotted 
in an RA-cos(<5) equal-area rectilinear projection. The angular 
mask for the quasars is both determined by the spectroscopic 
coverage, as well as the overlap with the photometric LRG sample. 




M, = -2.66 logio(Lboi) +79.36 , 



(3) 



where Lboi is in Watts, and we co nvert from the bj to i 
band using the empirical relations in [Richards etHI l|2006l ). 
For the LSTAR sample at the median redshift of 0.43, with 
M*,i(z = 0.43) — —23.04, this implies a bolometric luminos- 
ity Lboi > 10^*'^ W. Assuming that the quasars are radiating 
at the Eddington rate (I/Edd = lO^^'^ (Mbh/lO^^Mo) W), this 
implies black hole masses Mbh > 3 x 1O^M0. 

Table [1] summarizes the various subsamples we consider 
in this paper (discussed further in Sec. 13. 2p . 



Figure 3. The conditional distribution of absolute magnitude 
with redshift for our sample of QSOs. The lines plot the 16, 50 
and 84 per cent contours. The lower panel plots the absolute 
magnitude relative to M*, estimated from the 2dF and 2SLAQ 
survey. The flattening/upturn at low redshift is due to the Mi = 
—22.0 cut on the sample, to minimize contamination from the 
host galaxy. The vertical lines (as in Fig. [TJ show the redshift 
boundaries of our samples. 
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2.2 Luminous red galaxies 

We cross correlate the above quasars with a sample of 
luminous red galaxies (LRGs) selected from the SDSS 
imaging data. The sample selection, angular mask, and 
redshift distributions hav e been described in detail in 
IPadmanabhan et al.l l)2007i ). and we refer the reader to the 
details there. These galaxies have well characterized pho- 
tometric redshifts and errors {5z ~ 0.03), allowing us to 
deconvolve the photo metric redshift distribution to obtain 
the underlying dN /dz jPadmanabhan et al.ll2005l ). We con- 
sider LRGs with 0.25 < Zphoto < 0.55, trimmed to an angu- 
lar mask that covers 3528.04 deg.^ in the northern Galactic 
hemisphere; this results in a sample of 454882 LRGs. We 
then divide the sample into six redshift ranges, summarized 
in Table [21 of photometric redshift width Az = 0.05, where 
the redshift boundaries are chosen to select approximately 
homogeneous samples as the 4000 A break shifts from the 
g to the r-band. We found that the LRG properties varied 
significantly over the redshift range making fine Zphoto bins 
essential for proper modeling, the choice of Az = 0.05 be- 
ing determined by the photometric redshift errors. The LRG 
samples are summarized in Table [2j while the deconvolved 
redshift distributions are in plotted in Fig. [T] 



3 CLUSTERING 
3.1 LRG Clustering 

3.1.1 Methods 

We measure the clustering of the LRG sample in each of the 
six photometric redshif t slices using the angular clustering 
estimator described in IPadmanabhan. White fc EisensteinI 
l|2007l ). We define 

uj{es) = 2-K [ ede G{e,es)w{e) , (4) 

where 

x>l, (5) 

with X = 6/6s. As was shown in 

IPadmanabhan. White fc EisensteinI (|2007l ). this esti- 
mator partially deprojects the angular correlation function, 
yielding a robust estimate of the 3D real space correlation 
function on scales of ~ |x^s where x is the mean comoving 
distance to the redshift slice under consideration. 

We implement the above estimator using DD / RR—1 as 
our estimate of the a ngula r correlation function. Although 
the iLandv fc Szalavl (|l993l ) estimate is a more traditional 
choice, the contiguous wide-area coverage of the SDSS imag- 
ing makes it unnecessary in our case, and we choose the 
simpler estimator. Substituting this into Eq. |31 we obtain, 



= 2tt 



edo G{e,e,)^ , 



(6) 



where we use the fact that the area-weighted integral of 
G{x) vanishes by construction. In order to proceed, we note 
that the RR term is a purely geometric term determined by 
the survey mask. Although this is traditionally estimated by 
measuring random-random pairs in the same binning as the 



DD pairs, IPadmanabhan. White fc EisensteinI |2003) point 
out that on scales much smaller than the size of the survey, 
RR is described by 



cx 27r6'A6l $(61) 



(7) 



where '^{9) is a smooth function; we obtain a good fit to $ 
using a fifth-order polynomial. Having fit "I>, we can make 
our 9 bins arbitrarily small without incurring any Poisson 
noise penalty. This allows us to rewrite Eq. |S] as a weighted 
sum over pairs. 



oj{9s 



G{9„9s) 



iG-D-D 



<d{9s - 9) , 



(8) 



where O is the Heaviside step function. 

We estimate the covarian ce matrix of our meas urements 
by bootstrap resampling (e.g. lEfron fc Gonelll983l ). An im- 
portant advantage of uj{9s) is its insensitivity to clustering 
on scales ^29 s. This allows us to subdivide the survey into 
41 approximately filled spherical rectangles, 7.5° in the RA 
direction and 0.15 in the sin(5) direction. Computing uj{9s) 
for each of these subsamples yields 41 independent, iden- 
tically distributed realizations of lj{9s). The independence 
of the subsamples is a direct consequence of the estimator; 
this is not true for the more traditional Vli{9). We then esti- 
mate both the average and covariance matrix by bootstrap 
resampling these 41 realizations. In order to improve the nu- 
merical stability of this procedure, we scale why 9^, thereby 
removing the artificially large condition number of the co- 
variance matrix that arises due to the large dynamic range 
of uj. The resulting covariance matrix is very well-behaved, 
with no anomalously small modes that need be removed. 
Note that this is not the case for w(S) which is sensitive 
to large scale fluctuations, which in turn lead to unphysical 
modes in the covariance matrix that must be further con- 
ditioned. The insensitivity of uj{9s) to long wavelengths is 
an important advantage when attempting to estimate the 
covariance matrix from the data itself. 



3.1.2 Results 

The angular clustering of the LRGs in the six redshift slices 
is shown in Fig.|4l using the estimator described above; also 
plotted is the approximate physical scale probed by tj at a 
given 9s. Note that the errors between different points are 
correlated, and we use the full covariance matrix in all fits. 

Given the observed angular clustering, one can infer 
the underlying three dimensional clustering of the sam- 
ple; we do this using two methods. The first assumes that 
the LRG clustering traces the dark matter with a scale- 
indepe ndent bias on large scales. We compute uj{9s) assum- 
ing the lSmith et al.l l|2003l ) prescription for the shape of the 
nonlinear dark matter power spectrum, and estimate the 
large-scale bias by fitting to the data on scales larger than 
0.2°, corresponding to physical scales ;^2 /i~^Mpc. The best 
fit models are plotted i n Fig. [31 with the c orresponding bias 
values in Table [21 The [Smith et al.l (|2003l ) prescriptions de- 
viate from the observed clustering on small scales; as one 
might expect: LRGs do not trace the dark matter on these 
scales. 

The second method attempts to model the observed 
clustering by populating halos in a dark matter simulation 
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"LRG 


^photo-range 




-^modc 


oz 


n 




h 

'.'const 


"halo 


LRGl 


29,660 


0.25 <z < 0.30 


0.276 


0.287 


0.029 


0.87 


4.7 


1.75 ±0.05 


1.71 ±0.05 


LRG2 


32,527 


0.30 < z < 0.35 


0.326 


0.312 


0.033 


0.86 


4.1 


1.77 ±0.06 


1.77 ±0.05 


LRG3 


41,051 


0.35 < z < 0.40 


0.376 


0.375 


0.049 


0.84 


3.4 


2.36 ±0.05 


2.15 ±0.07 


LRG4 


60,294 


0.40 < z < 0.45 


0.445 


0.452 


0.058 


0.81 


3.7 


2.28 ±0.05 


2.09 ±0.05 


LRG5 


104,131 


0.45 < z < 0.50 


0.506 


0.488 


0.048 


0.79 


4.7 


2.02 ±0.04 


1.90 ±0.04 


LRG6 


95,605 


0.50 < z < 0.55 


0.552 


0.541 


0.051 


0.78 


4.2 


1.90 ±0.05 


1.76 ±0.05 



Table 2. Properties of the LRG samples, showing the number, photo-z range, mean and modal redshift, width, growth factor (normalized 
to unity today) at the modal redshift, number density (in 10^* /I'^Mpc"^) and large-scale bias. We estimate the large-scale bias by fitting 
to w(5s) assuming scale-independent bias or fitting to a halo model (see text). The quoted errors are purely statistical. The biases 
estimated from the different methods have systematic errors at the 5 per cent level. We adopt fehalo as our fiducial value. 



with LRGs; we refer the reader to Appendix [SI for details. 
The best fit a; are in Fig. U] and the predicted large scale 
bias values are in Table [2] Taking into account the likelihood 
of scale-dependent bias and the 5 per cent systematic uncer- 
tainty in the halo modeling (see Appendix [A} these are in 
reasonable agreement with the simple fits described above. 
The halo models (Fig. [5} reproduce the observed correla- 
tions well, including the prominent break in the correlation 
function (the best fit x'^ values are in Table . Our LRG 
samples populate a broad range of halo masses, with an ap- 
proximate power law dependence of the mean number of 
LRGs per halo with the halo mass. Furthermore, we find 
that halos with masses ~ W^Mq have one LRG in them. 
An important by-product of this process is that we obtain 
mock realizations of the LRG sample; we use these below to 
interpret the clustering measurements of the quasars. 

The LRG clustering amplitude is consistent with being 
non-evolving with redshift, implying a bias that evolves as 
h(z) ~ 1/D[z). An exception are LRG3 and LRG4, which 
have a significantly higher bias. The reasons for this were dis- 
cussed in detail by IPadmanabhan et al.l (|2007l ). Comparing 
the bias values and HOD fits, we see that the LRGs can be 
conveniently grouped into three slices of width dzphoto = 0.1, 
each of which samples a homogeneous population of galax- 
ies. 

Finally, we note our re sults are consistent with those of 
IPadmanabhan et al.l |20o3) correcting for differences in the 
fiducial cosmologies. That paper also performed a number 
of systematic tests on the LRG sample, and we simply refer 
the reader to that work instead of repeating them here. 

3.2 QSO Clustering 

3.2.1 Methods 

The quasars in our sample have spectroscopic redshifts and 
thus we know (up to small uncertainties due to peculiar 
velocities) a physical distance to each object. This allows 
us to work in terms of transverse separation rather than 
angular separation, i.e. to measure 



Wp{R) 



(9) 



rather than 'w{9) or its generalizations. As with the LRGs, 
this is a real space measurement, avoiding the need to model 
redshift space distortions. 

We start by considering quasars in a narrow redshift 
range (a comoving distance xo away), and correlate them 



with LRGs with a normalized radial distribution, /(x)- In 
the fiat sky approximation, the angular correlation function 
is given by. 



w{e) 



dx fix) ^ (VxlO' + ix-Xor) 



(10) 



where the usual second integral over the quasar redshift 
distribution has been eliminated because the quas ars have 
spectroscopic redshifts. We can now make the usual [Limber! 
(1 953 ) approximation for xod much smaller than the scales 
over which f{x) varies. This allows us to hold fix) fixed at 
fixo) in the integral, yielding 

wie) ~ fixo) J dx C (Vxge^ + (x-xo)^) (11) 

= /(Xo)^p(i?) . (12) 

Note the second use of the assumption of a peaked integral 
to re-extend the limits of the integral to ±oo. The second 
equality just recognizes the integral as the projected corre- 
lation function Wp at transverse separation R — xoS- 

To generalize to a broad redshift slice, we start with 
the standard assumption that Wp{R) does not evolve over 
the slice. One then has two choices - the first is to esti- 
mate 1^(9) over narrow redshift slices, estimate Wp(i?) for 
each of these slices using Eq. [T2]and then average. This has 
the disadvantage that each individual WpiR) measurement 
is extremely noisy, and potentially sensitive to noise in fix)- 
The other approach (which we adopt) makes use of the fact 
that 9 in Eq. [T^] is simply a label, and can just as easily 
be replaced by R. This has a simple interpretation - one 
computes the angular correlation function over the broad 
redshift range but replaces the angular separation by the 
transverse separation R computed assuming that the LRGs 
are the same redshift as the quasar being correlated with. 
We use the estimator, 



weiR) 



QGjR) 
QRiR) 



(13) 



where QG and QR are the quasar-galaxy and quasar- 
random pairs, and our notation makes explicit that we are 
binning in physical transverse separation. Since we assume 
that WpiR) is constant over the redshift range. Eg. 1121 yields 



MR) = {fix))MR) , 



(14) 



where the average is done over the quasar redshift distribu- 
tion. Note that this formulation addresses both the prob- 
lems of the first implementation. Finally, we point out that. 
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X eJZ (Mpc/h) X e./2 (Mpc/h) X (Mpc/h) 




0.1 1 0.1 1 0.1 1 



e. (deg.) e, (deg.) e, (deg.) 



X eJZ (Mpc/h) x 9./2 (Mpc/h) X e./2 (Mpc/h) 

1 10 1 10 I 10^ 




0.1 1 0.1 1 0.1 1 



e, (deg.) e, (deg.) 9, (deg.) 



Figure 4. The observed u) for the six LRG redshift shoes, as a function of the filter scale 9s ■ Recall that uj probes the angular correlation 
functi on on scales ~ ds/2; the corresponding physical scales are also shown. We also plot the best fit models, both from the Smith ct all 
l|2003l) fitting formula [dashed/red] for the nonlinear dark matter clustering, as well as from our halo model fits [solid/blue]. The dotted 
vertical line marks the angular scale beyond which we fit the HaloFit models. These fits deviate from the observed clustering on small 
scales; LRGs are not distributed like the dark matter on these scales. The halo model correlation functions are estimated from the same 
realizations in the 500/i~^Mpc simulation that we use to interpret the quasar-LRG cross-correlations. Note that the data are well fit by 
the halo model on both small and large scales. 



in Eq. 1131 one only requires the angular selection function 
of the LRGs to estimate QR. 

Note that one could use the estimator described by 
Eq. |3]to measure the cross-correlations. However, the lower 
signal-to-noise of the quasar-LRG cross correlations elim- 
inates the advantages of the estimator, and therefore, we 
choose the simpler estimator. 

Estimating the covariance matrix for our sample is sim- 
plified by the low areal density of the quasars, making them 
effectively independent for the scales of interest. We there- 
fore estimate the covariance matrix by simply bootstrapping 
the individual quasars. As with the LRGs, we remove the 
artificially large condition number by scaling Wp{R) by R 
before estimating the covariance matrix. Finally, as a check, 
we note that we obtain consistent results if we replace the 
bootstrap covariance matrix with a jackknife estimate. 

3.2.2 Results 

In order to determine what subsamples to cross-correlate 
with, we start with the observation (see Sec. 13. ip that the 
LRGs can be grouped into three homogeneous slices of width 
0.1 in photometric redshift, i.e. grouping slices (1,2), (3,4), 
and (5,6) of Table [2] together. We use the mean redshift of 
each LRG redshift slice and its width to determine the red- 
shift range of the quasars to cross-correlate with; this defines 



the redshift ranges of the quasar subsamples in Table[T] Note 
that since we use the true redshift distribution of the LRGs 
(as opposed to the photometric redshift distribution) to de- 
termine the mean and width of the slices, we automatically 
correct for any biases and asymmetries in the photometric 
redshift errors. For each of these three quasar redshift slices 
(denoted as ALL below), we further consider the following 
subsamples - the LSTAR sample defined in Sec. [2j restrict- 
ing to quasars targeted with the latest version of the SDSS 
target selection algorithm, and a bright and faint subsample 
split at the median luminosity for each redshift slice. 

As anticipated earlier, we find that restricting to the lat- 
est version of the SDSS target selection algorithm gives the 
same cross-clustering power against the LRGs as the ALL 
sample. Since it has better statistics we will use the ALL 
sample below. Furthermore, as is evident from Fig. [S] the 
luminosity baseline is rather small, and no clear trend with 
luminosity emerges, given our errors. We therefore present 
detailed results for only the ALL and LSTAR samples below. 
The higher redshift slices do not extend significantly below 
L^, so to simplify the interpretation we use the LSTAR sam- 
ple as our fiducial sample. Fig.|6]plots the cross-correlations 
for the LSTAR sample. 

We present both power law and large scale bias fits to 
both these subsamples. In order to estimate the mean and 
error for any parameter p (possibly a vector), we use 
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Figure 5. Best fit HODs for the six LRG slices we consider in this paper. The shaded region denotes the errors, as estimated by Monte 
Carlo. 



(P> 



and 



/ dpC{p)p 



(Pi) - {P^) 



(15) 



(16) 



where the likelihood is defined hy C = exp (— with 
computed using the full covariance matrix. For the power 
law fits, we adopt a two parameter model, 



Wp{R) 
R 



/;Fr[(7- l)/2] /roV 



r[7/2] 



m 



(17) 



which corresponds to a 3D cross-correlation of the form 
^(r) = (r/ro) ' . We fit this model to the measured cor- 
relations on all scales. In order to determine the large-scale 
bias we compute ■Wp{R) for the dark matter at z = 0.3 , 
0.42 and 0.53 using the prescription in ISmith et al. I (|2003D . 
This is then scaled by an ii-independent multiplier to ob- 
tain the best fit to the data in the range R > 2 h~^Mpc. In 
the limit of scale-independent, deterministic bias the multi- 
plier is bgbLRG (see Table [3] for the values). We transform 
these into bg by using bLRG for slices 1, 3 and 5, scaled to 
the corresponding redshift by the growth factor. Note that 
the LRG clustering amplitude is close to constant with red- 
shift, so we would obtain consistent numbers if we had used 



the other slices. The lower limit in the fit, 2/i~^Mpc, was 
determined by the scale at which Wp{R) from QSO-LRG 
cross-correlations in our mock catalogs (see ii4.2|) showed 
significant scale-dependent bias. 

The results for the large scale bias and power law fits 
are in Table [31 while Fig. [7] plots the evolution of the clus- 
tering amplitude of the LSTAR sample as a function of red- 
shift. Our results are consistent with a constant clustering 
amplitude from z = 0.25 to 0.6, corresponding to a bias of 
1.09 ± 0.15 at z = 0.43. 

Table [4] summarizes our results compared with previ- 
ously published workQ Our results strongly favour the gen- 
eral consensus that the bias of low redshift quasars is ^ 1 ; 
this is also consistent with the models of iHopkins et al 



(2007) as well as the previous extrapolations bv lCroom et al 



Mvers et al 



(2005j). There are two significant exceptions 
(2007a|) find 1.93 ± 0.14 based on a photometrically se- 
lected sample of quasars. It is possible that contamina- 
tion by a high redshift population could boost the mea- 
sured bias values. The more intriguing discrepancy is with 



^ We do caution the reader that the errors for a number of these 
measurements are simply Poisson errors, and ignore correlations 
between different scales, and are therefore likely underestimated. 
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iMountrichas et al] l|2008l ) who analyze a similar sample to 
ours, also in cross-correlation with LRGs, and find biases be- 
tween 1.90±0.16 and 1.45±0.11 depending on the particular 
LRG and quasar sample they cross-cor relate again s t. Th ese 
results are also discrepant with ida Angela et alj (|2008| ) - 
with whom we are consistent - who analyze the same sample 
in auto-correlations. Furth ermore, the scatter in the differ- 
ent subsamples analyzed bv lMountrichas et al.l (|2008h signif- 
icantly exceeds their quoted errors, suggesting either a sys- 
tematic in their analysis, or an underestimate of their errors. 
Using the observed scatter between the different subsamples 
as an estimate of the error yields a value consistent with 
our mea surement. Finally, t wo r esults not presented in Ta- 
ble Ha re [Serber et al.1 (120061 ) and lStrand. Brunner fc MversI 
l|2007fl . There are a number of subtleties with interpreting 
these results (resulting in misunderstandings in the litera- 
ture) ; we therefore defer a detailed discussion of these results 
to Appendix [B] 



4 INTERPRETATIONS 
4.1 The Large Scale Bias 

The large-scale bias of any population of objects provides 
information on the mean dark matter halos mass hosting 
that population. Specifically, if N{Mh) is the mean number 
of QSOs hosted by a halo of mass Mh then 



(b) 



alvih 



(18) 
(19) 



where duh/dMh is the (com oving) number density of ha- 
los per mass interval (e.g. ISheth fc TormenI Il999l ) and 
bh(Mh) is the bias a ss ociated with halos of that mass (e.g. 
Cole fc Kaiseil Il989l : ISheth fc TormenI 1 19991 : ISheth et al.l 
200 ll ). Note that (b) is independent of the normalization 
of N{Mh)- Also, since observable quasar activity is a tran- 
sient property, the observed quasar number density fiqso 
depends upon the average duty cycle, /on, 



nqso = ,fonn 



(20) 



To set the scale. Fig. [7] compares the observed large- 
scale bias with that of halos of fixed mass in our as- 
sumed cosmology. Our data are consistent with host ha- 
los having a ma ss 10^^'^ — lO^'^'^ h~^Mp), in agreement 
with earlier work (IPorciani. Magliocchetti fc Norbere 20041: 
Croqm et all l2005l : iPorciani fc Norberd l2006l : iLidz etal l 
20061 ). PI Note that our constraints are significantly stronger 
for higher as opposed to lower halo masses, since bh{M) is 
a rapidly rising function above ~ 10^^ h~^MQ but slowly 
asymptotes to a constant (6 ~ 0.5) for lower masses. As we 
see no evidence for evolution in the QSO bias, we use the 
mean value from all three slices below. 

Ifowever we don't expect QSOs to inhabit halos of a 
single mass. To place constraints on the range of halos in 
which quasars may be active we consider two illustrative 

^ A little appreciated uncertainty in this conversion comes from 
differences in fitting functions to (M) resulting in an additional 
error of 50 per cent (0.2 dex) in mass. 



models. First we imagine that QSOs brighter than live 
in halos more massive than Mmin, with each halo above 
Mmin hosting exactly one QSO with probability /on. This 
gives N{Mh) = Q{Mh - Mmin), where the Heaviside 6 
function is unity for positive arguments and zero otherwise. 
Fig.[S]and Table [S] summarize the constraints on this model. 
As anticipated by our simple scaling argument above, our 
measurements suggest Mmin ~ 10^^ h~^MQ corresponding 
to an average halo mass (M) ~ lO^^/i-^M0. We further 
strongly disfavour models with Mmin S> W^^ Mq. Our 
lower mass limits are significantly weaker for the reasons 
discussed above, requiring Afmin > 10'' '' h-^MQ, with even 
lower masses still providing marginally acceptable fits. How- 
ever, such masse s are disfavoured by the l ocally observed 
Mbij — o relation l|Ferrarese fc Merrittll2000l : iTremaine et al.l 
2002). Assuming an Eddington limited accretion rate, the 
LSTAR sample should be powered by black holes with 
masses M;,^ > IQ^ Mq, which live in bulges with a ~ 
lOOkm/s or Mbuigc a few times 10^ h~^MQ. Given mass-to- 
light ratios of a few, this disfavours Mmin < 10^° h~^MQ. 

The second model we consider starts from the assump- 
tion that QSOs cluster like a random sampling of a luminos- 
ity or color subsample of galaxies. This motivates a form. 



N{Mh) = e(Af - Afmin) 



1 + 



V20A/mi„y 



(21) 



found to be a good description of galaxies at both low and 
high redshifts (e.g . IZehavi et all l2005l : IConrov et all l2006l : 
IWhite et al.l 12007 ) . This has a central galaxy in all halos 
above Mmin and on average (M/20Mmin)°' satellites in each 
halo, where the factor of 20 is inspired by fits to SDSS galax- 
ies and N-body simulations. We consider cases with a = 0.5 
and 0.9; the former to model blue galaxies (which are under- 
respresented in high-mass halos) and the latter to describe a 
luminosity selected sample. (The models are supposed to be 
illustrative.) As before, a fraction /on of the QSOs are "on" 
at the time of the observation. Fig. [8| and Table [H again 
summarize the constraints from the large-scale bias, which 
are very similar to those obtained for the threshold models. 
Within the context of this model we can translate this into 
a constraint on the space density of quasar hosts and hence 
on their luminosity. The preferred value of Mmin suggests 
n ~ 5 X 10~^ M pc~^ or a lurninosi ty L < 0.1 using the 
blue-galaxy LF of IPaber et al] (|2007l '). The 95 per cent con- 
fidence level upper limit on b gives n ~ 4 x 10~^ Mpc~'^ 
or a luminosity L ^ L^,. This suggests these low redshift 
quasars live in relatively faint galaxies. The number density 
at fixed bias is lowered if we allow scatter in the L — M 
relation (i.e. a smooth turn on in Eg. I21|l . as is likely. With 
log- normal scatter in M at fixed L of uinM = 1 the space 
density is reduced a factor of ~ 4 at fixed bias and the up- 
per limit on the threshold lu minosity becomes 1.1 (c.f. 
I White . Martini fc Cohn"2008'). Using the stellar mass func- 
tions in Bundy ct al. (2007), we find that these number den- 
sities correspond to stellar masses M* < 2 x lG^^h~^ Mq. Of 
course QSO hosts may not inhabit halos in the manner as- 
sumed by Eq. 1211 If QSO hosts are under-represented in 
intermediate mass halos then it is possible to have a lower 
number density and fo ~ 1. This is what is seen in e.g. a sam- 
ple of galaxies with Mb < —21 and a star formation rate 
> IMpyr"^ in the Millennium simulation (ISpringel et al.l 
I2OO5I : ICroton et all I2OO6I : iDe Lucia fc Blaizotl l2007l ). The 
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Sample 


LRG Sample 


bQSobhRG 


bqso 




ro (Mpc/h) 


7 


(ro) (Mpc/h) 


(7> 


r 




ALLO 


1 


2.73 ±0.45 


1.60 ±0.26 


1.78 


6.70 


1.71 


6.48 ± 0.63 


1.70 ± 0.10 


0.18 


8.01 


ALLl 


3 


2.56 ± 0.47 


1.17 ±0.22 


4.03 


6.44 


1.82 


6.29 ± 0.54 


1.81 ± 0.09 


-0.30 


6.70 


ALL2 


5 


2.30 ±0.38 


1.20 ±0.19 


2.86 


5.56 


1.82 


5.46 ± 0.42 


1.82 ± 0.08 


-0.28 


5.94 


LSTARO 


1 


2.50 ±0.62 


1.46 ±0.37 


3.88 


6.60 


1.83 


6.23 ± 0.90 


1.82 ± 0.15 


-0.13 


7.42 


LSTARl 


3 


2.12 ± 0.64 


0.96 ±0.29 


4.54 


6.00 


1.93 


5.75 ± 0.79 


1.92 ± 0.14 


-0.49 


10.22 


LSTAR2 


5 


2.08 ± 0.41 


1.08 ±0.20 


5.15 


5.42 


1.91 


5.29 ± 0.48 


1.91 ± 0.09 


-0.44 


8.68 



Table 3. The large-scale quasar bias and power law fits for the ALL and LSTAR samples. We fit scales _R > 2h~^Mpc (5 points) 
for the bias, and all scales (12 points) for the power law model. The second column lists the LRG slice assumed for the bias (see text 
and Table |2] for details) , while the third and fourth columns list t he amplitude of th e quasar-LRG cross-correlation and the implied 
(scale-independent) bias, assuming a dark matter ^(r) given by the ISmith et aP l|2003l ) model. Note that we ignore the contribution of 
the error in the LRG bias in the derived QSO bias. The best fit ro and 7 values are listed under ro and 7, while likelihood-averaged 
values are under (ro) and (7). Note that the errors in ro and 7 are correlated; the cross-correlation coefficient is under r. 




Figure 6. The observed Wp{R) measuring the cross-correlation between quasars and the LRG samples. The solid [red] line shows the 
HaloFit prescription for the nonlinear dark matter correlation function, normalized by a scale independent bias to best fit the observed 
correlations. The dashed [blue] line shows the best fit power law model. 



z 


-^min 


bqso 


Reference 


0.25 < 2 < 0.6 


L* 


1.09 ±0.15 


(1) 


0.3 < 2 < 0.68 


0.4L* 


1.27 ±0.20 


(2) 


0.4 < 2 < 1.0 


O.IL* 


1.93 ±0.14 


(3) 


0.7 < 2 < 1.4 


O.IL* 


1.09 ±0.29 


(4) 


z ~ 0.6 


0.4L* 


1.10 ±0.20 


(5) 


z ~ 0.6 


0.4L* 


1.90 ±0.16 


(6) 


z ~ 0.6 


2.5L* 


1.45 ±0.11 


(7) 


z < 0.3 


0.4L* 


0.97 ±0.05 


(8) 



Table 4. A summary of previous low redshift quasar clustering 
results, compared with results in this work, scaled to the cos- 
mology assumed her e. (1) This work (2) | Croom et all (|2005|) (3) 
iMvers et all ll2007ah (4) ICoil et al.l ||2007|) (5) Ida Angela et al.l 
||2008| ') (6.7) iMountrichas et al.l 1I2OO8I ') (8) ICroom et al.l 1120011 
For the results from lColTet'^ 11200^, we scale the relative bias 
presente d there by the lar ge scale bias b = 1.22 of all DEEP2 
galaxies l lZheng et al.ll2007^ 



HOD of these galaxies is approximately described by a log- 
normal distribution peaking at M ~ 10^^ AIq, with a 
power law distribution at high masses. There is however a 
deficit of galaxies at ~ 10^^ AIq; this allows one to have 
a low number density, without a high bias. Note that while 
our data are unable to constrain such flexible models, a large 
fraction of this uncertainty derives from the fact that the 



HOD 




(M) 


A/95 
mm 


Thresh. 


0.3 


1.9 


2.0 


a = 0.5 


0.3 


3.1 


1.6 


a = 0.9 


0.1 


6.1 


1.3 



Table 5. Halo masses derived from the average large-scale bias of 
our QSO sample. The masses are quoted in units of lO^'^ h~^MQ. 
The different HOD models are described in detail in the text. The 
superscript 95 refers to the 95 per cent c.l. upper limits. 

QSO number density does not add any constraints to the 
HOD. 

Given the above models and caveats, and the ob- 
served space density of quasars, we constrain the duty cy- 
cle fEq. | 20|) to be < 0(10"^), consistent with estimates by 
iDunlop et all (|2003l ) from the luminosity function. Convert- 
ing the duty cycle, /on, into a lifetime is somewhat ill de- 
fined. If we assume Iq = fontn, with tn the Hubble time, we 
find tg < 10^ yr. These lifetimes are broadly consistent with 
those derived at 2 ~ 2. On the other hand, the Hubble time 
is significantly longer, and the duty cycles are significantly 
lower. 

The exact timescale to use in the above conversion is 
not well defined; our choice of the Hu bble time for the halo 
lifetime is an approximation (see e.g. iMartini fc Weinberg! 
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Figure 7. The bias of the three LSTAR quasar subsamplos 
[squares] as a function of redshift, as well as the mean bias [circle]. 
We consider the bias scaled by the growth factor to focus on the 
evolution of the clustering amplitude. Also plotted [dashed lines] 
is the bias of halos with masses (from bottom to top) ranging from 
logio(M/M0) of 10.5 to 13.5 in steps of dlogio(M/MQ) = 0.5. 
Note that the mean large scale bias of 1.09 ± 0.15 at z = 0.43 
suggests that quasars live in halos of mass ^ W^'^h~^ Mq. 




Figure 8. Afmin for the threshold [solid, black] and galaxy-like 
a = 0.5 [dotted, red] and a = 0.9 [dashed, blue] HODs, as con- 
strained by the average large scale bias of the LSTAR sample. The 
dotted hatched region marks the region disfavoured by quasar 
lifetimes and the Mj^j^ — a relation. The dashed hatched region is 
excluded at 95 per cent confidence. See the text for more details. 



Figure 9. Comparison of the cross-correlation of simulated 
LSTAR2 and LRG5 samples with the observations. The solid 
(x^ = 12.7) and dashed (x^ = 8.1) lines assume a galaxy-type 
N{Mi^) with a = 0.3 and 0.9 respectively, while the dotted 
lines (bottom to top) assume a threshold N(M}^) with a mass- 
independent satellite fraction of 0,0.5 and 1 (x^ = 19.8, 15.0 and 
11.7 respectively). We fix M^^^ to the best fit values in Table [S] 
note that the large scale clustering for all the models are identi- 
cal. The x^ values are computed using all 12 points and the full 
covariance matrix. 

l200ll . for a different approximation) . The differences between 
tfie various choices are of tlie same order of magnitude as 
tfie systematics in modeling the quasar host number density, 
and we tlierefore opt for simplicity. We however caution the 
reader that these numbers should only be treated as order- 
of-magnitude estimates. 

If QSOs are radiating at the Eddington limit Ledd then 
the minimum A^bh in our LSTAR sample is 3 x lO^'M©. This 
value is consistent with the estimates from the Mbh — Mhaio 
relation (Fcrrarcsc'2002'), Mbh ~ 2 x lO'' - 3 x IQ^Mq, with 
the differences coming from different assumptions about the 
halo profiles. This suggests L/Ledd ^ 0.0 1 — 1, consistent 
with the results of ICroorn et all (|2005h and Ida Angela et al.l 
(2008), although we find no evidence of super-Eddington 
accretion. 



4.2 Small Scale Clustering 

Due to the larger number density of LRGs with which we 
have cross-correlated our QSOs we are able to measure 
the clustering down to smaller scales than would otherwise 
have been possible. Interpreting the small-scale clustering 
is however complicated by the uncertain relation between 
the galaxies which host active QSOs, the galaxies which 
are selected as LRGs and their parent host halos. A full 
interpretation would require knowledge of the joint distri- 
bution P{Nqso, NlrgIMii), which cannot be meaningfully 
constrained with the limited data we have available. It is 
however straightforward, with the aid of the mock catalogs 
described previously, to predict the cross-clustering for any 
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well specified QSO scenario. In this section we examine the 
illustrative models introduced in the last section, focussing 
on the z ~ 0.5 data for definiteness. 

We begin with the model in which QSOs are hosted ex- 
clusively by the central galaxies of halos above some thresh- 
old Mmin. As with all of the models we shall ignore the 
number density constraint by postulating that a random 
fraction, /on, of the possible hosts are seen as QSOs at any 
given time. If the probability of a host being on is indepen- 
dent of the host properties the clustering is unchanged and 
we can use the host population - with its better statistics 
- to compute Wp. We find that this model cannot simulta- 
neously account for the measured large-scale bias and the 
large amplitude of the small-scale correlation function (see 
the lowest dotted line in Fig. |9} . To fit the former the mean 
QSO-weighted halo mass must be low, but such low mass 
halos contain only central LRGs, not satellites, so there are 
no QSO-LRG pairs with separations 0(100 h~^kpc). We can 
keep the halo occupancy the same but distribute the QSOs 
within the halos as both centrals and satellites. We assume 
satellites follow an NFW profile l|Navarro. Frenk &: Whit3 
1 19961 ) with a concentration of 10. - the precise details don't 
matter for our purposes, other choices produce qualitatively 
similar results and trends. By making 50 per cent of the 
QSOs satellites we boost the power on small scales (Fig. |9]), 
as it is now possible to have QSO-LRG pairs in the smaller, 
lower mass halos. 

We can also take the form of Eq. 1211 imagining that a 
random sample of galaxies will host QSOs at the time of 
observation. We place a central galaxy in all halos above 
Mmin and place a Poisson number of satellites, with mean 
(M/20Mn,in)", distributed like a c = 10 NFW profile. The 
large-scale bias is again set by the mean halo mass, which 
is larger for larger a at fixed Mmin- As a is increased the 
spread in small-scale clustering amplitude with Mmin de- 
creases, with models lying very close to the data. A model 
with a ~ 0.9 (solid line in Fig. (9)1 gives very good fits, with 
X'^/dof < 1, when Mmin ^ 10^^ Mq but lower values 
of a (dashed line in Fig. ^ are not excluded. Due to the 
strong covariance between the Wp{R) points, the constraint 
on Mmin from the full model is not stronger than that from 
just the large scale points. 



5 CONCLUSIONS 

We measure the small scale clustering of a sample of 
~ 400, 000 photometric luminous red galaxies and their 
clustering around a volume limited sample of ~ 2000 z < 0.6 
low redshift QSOs. By using a new statistical estimator, we 
are able to obtain precise measurements of the LRG angular 
correlation function, which coupled with their precise and 
well-characterized photometric redshifts, allowed us to con- 
strain how LRGs populate dark matter halos. We find that 
LRGs have a clustering amplitude that is consistent with not 
evolving with redshift, and corresponding to a large scale 
bias b 2 at z — 0.5. The best fit halo occupation models 
suggest that they occupy halos > 10^^ Mq, with approx- 
imately one LRG in every 10^^ Mq halo. We use these 
halo occupation distributions to construct mock catalogs of 
LRGs. Attempting to match the observed cross correlation 
of LRGs with QSOs by populating these same mock catalogs 



with QSOs allowed us, for the first time, to start to probe 
how quasars inhabit dark matter halos. 

The cross-correlation of QSOs and LRGs is well de- 
scribed on all measured scales by a power law of slope 
~ 1.8 ±0.1 and a scale length of ~ 6 ± 0.5 /i~^Mpc, consis- 
tent with observed slopes and amplitudes for local galaxies. 
It is also well described by the nonlinear matter correla- 
tion function, scaled by a constant bias, although there is 
some evidence for deviations from this form at the smallest 
scales. Such a deviation is however not unexpected, and is 
seen for most galaxy samples, which are better described 
by power laws down to small scales. S ince this i s in appar- 
ent contradiction with the results of Se rber et al.l (2006) and 
[strand. Brunner fc MversI 120071), we revisit their measure- 
ment within the framework of cross-correlations developed 
in this paper (Appendix [Bj . We explicitly show that their 
data are fit by a power law of slope 1.9, with no deviations 
on small scales, and that the claims of an excess come from 
subtleties in interpreting their measurements. 

The large scale bias b — 1.09 ± 0.15 is consistent with 
most previous measurement s and th eoret ical models, the ex- 
ccptions being iMvers et al.l |2007a) and iMountrichas et al] 
(2008)i possible reasons for this discrepancy are discussed in 
ij3.2.2l We see no evidence for variations of the bias with red- 
shift or luminosity. The observed large scale bias constrains 
quasars to reside in halos with a mean mass of 10^^ Mq, 
with uncertainties of a factor of a few from the details of 
how the halos are actually populated. Our constraints on 
the halo mass are significantly stronger from above than be- 
low, since these halos are below the nonlinear mass scale 
and occupy the slowly varying region of the halo bias curve. 
This should be contrasted with measurements at higher red- 
shifts; even though the characteristic halo mass is the same 
at these redshifts l|da Angela et al.]|2008l ). it is now higher 
than the nonlinear mass and probes the steeply rising part 
of the halo bias curve. This problem is exacerbated when 
one considers realistic models of halo occupation. 

The mean halo mass can, in t urn, be used to con- 
strain the lifetimes of these QSOs (ICole fc Kaiser! Il989l : 
iHaiman fc HuH 1200 ll : iMartini fc Weinberg! |2001^ to be < 
10^ yr. This is consistent with measurements at high red- 
shift; on the other hand, the Hubble time is significantly 
longer today, and therefore the duty cycles are significantly 
shorter. We also discuss some of the theoretical uncertain- 
ties in estimating the number density of quasar hosts (and 
therefore the duty cycles and lifetimes), that arise from the 
difficulties in constraining how quasars populate dark mat- 
ter halos. We find the number densities of quasar hosts are 
only certain at the order of magnitude level, obviating any 
need for the detailed modeling (popular in the literature) of 
the conversion of number densities into qu asar lifetimes. 

Assuming local Mbii — Mhaio relations (|Ferrarese!l2002l ) , 
we estimate Eddington ratios between 0.01 to 1. We do not 
find any need for su p er-Ed dington accretion, in contrast 
with Ida Angela et all (|2008l ) who require super-Eddington 
ratios under certain assumptions for the halo profile. How- 
ever, the errors on both measurements are large, and are 
therefore consistent with each other. These Eddington ratios 
are also consistent with measurements at z ~ 2, suggesting 
no evolution in the Eddington ratio with redshift. 

Given our detailed modeling of the LRG clustering and 
the associated mock catalogs, we attempt to use the small 
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scale cross-correlation to constrain how quasars must pop- 
ulate dark matter halos. The size of our errors precludes 
being able to place interesting constraints on general quasar 
models. However, the forward modeling problem - taking a 
particular quasar model, and comparing it with our data - is 
straightforward. We find that our small scale measurements 
are inconsistent with quasars being a random subsample of 
all halo centers above a certain mass threshold, but become 
consistent if we assume > 25% of the quasars are satellites. 
We also find that our data are extremely well fit, if we as- 
sume that quasars are a random subsampling of luminosity- 
thresholded sample of galaxies, for luminosity thresholds 
between O.IL* to L*. The above results suggest that the 
host galaxies have a number density < 10~^/i"^Mpc~"^, corre- 
sponding to stellar masses < 10^^ Mq\ the exact val- 
ues are however sensitive to the particular choice of model. 

A second purpose of this paper was to demonstrate the 
modeling of how quasars populate dark matter halos. An in- 
teresting question therefore is how to optimize future mea- 
surements to gain the most leverage on this. Modeling QSO- 
galaxy cross-correlations differs from the traditional galaxy- 
galaxy auto-correlations in two important ways. The first is 
the lack of a constraint on the number density of the under- 
lying population that quasars are assumed to sample; which 
usually puts a strong constraint on the minimum halo mass. 
For quasars such a constraint must come from the bias and 
as discussed above this is only strongly constraining when 
one is on the steeply rising part of the halo bias curve. The 
second difference is that the one-halo term is only probed 
where the QSOs and galaxies occupy the same halos. For 
the quasar and LRG sample presented here, the mean halo 
mass for the quasars probes only the tail of the LRG HOD, 
making the constraints weaker than one might naively ex- 
pect. On the other hand, we emphasize that the accurate and 
well-characterized photometric redshifts of the LRGs were 
an essential prerequisite for doing the modeling in the first 
place - this could not have been done with the full SDSS 
photometric sample. This suggests that the ideal sample to 
correlate with would be a fainter sample of red galaxies. 
Furthermore, the modeling would be significantly easier at 
higher redshifts, where the quasar bias is higher. Various 
combinations of these will be available with the next gener- 
ation of imaging surveys, making it possible to significantly 
improve on the constraints presented here. 
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APPENDIX A: HALO MODEL 

In order to understand the manner in which luminous red 
galaxies and quasars inhabit dark matter halos wc make use 
of th e halo model (for a review see, e.g. ^Coorav & Shct^ 
120021 ). Within this formalism an accurate prediction of 
galaxy clustering requires a knowledge of the occupation 
distribution of objects in halos (the HOD) and their spatial 
distribution. In combination with ingredients from N-body 
simulations a specified HOD makes strong predictions about 
a wide array of galaxy clustering statistics. 

Our modeling of galaxy clustering is based on mock 
catalogs constructed within the HOD framework by popu- 
lating halos in a cosmological N-body simulation. We use a 
high resolution simulation of a ACDM cosmology {Q.m = 
0.25 = 1 ~ ^A, = 0.043, h = 0.72, = 0.97 and 
(Tg = 0.8). The linear theory power spectrum was computed 
by evolution of the coupled Einstei n, fluid and Boltzman n 
equations using the code described in lWhite fc Scot'O ll 19951). 
This code agr ees well with CMB fast ( Seliak fc Zaldarriagal 
Il99(|). see e.g. lSeliak et al.l l|2003l ). The simulation employed 
1024"^ particles of mass 8 x 10® /i~^Mq in a periodic cube 



of side 500/i"^Mpc using a TreePM code (White 2 0oi; for 
a com parison with other N-body codes see lHeitmann et al.l 
l2007h . The Plummer equivalent softening was 18/i"^kpc (co- 
moving). To check for finite volume and force resolution ef- 
fects we also looked at simulations of the same cosmology, 
with the same number of particles, in boxes 250/i"^Mpc and 
l/i~'^Gpc. The 250/i~^Mpc box turned out to be too small 
to model LRG clustering. 

For each output we generate a c atalog of halos usin g the 
Friends-of-Friends (FoF) algorithm (iDavis et al.l ll985') with 
a linking length of 0.168 times the mean inter-particle spac- 
ing. This procedure partitions the particles into equivalence 
classes by linking together all particles separated by less 
than a distance b, with a density of roughly p > 3/(27rfe'') ~ 
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100 times the background density. A comparison of the mass 
functions in the 500/i~^Mpc and l/i~^Gpc boxes suggests 
that at low particle numbers FoF tends to over-count the 
number of halos. To make the mass functions match in the 
overlap region we adjusted the mass of the halo s downward 
by a factor l — n~!^^^ (for a similar correction see lLukic et al.l 
I2OO7I ). 

To make mock catalogs we use a halo model which dis- 
tinguishes between central and satellite galaxies. We choose 
a mean occupancy of halos: N{M) = {Afgai(A'/haio)). Each 
halo either hosts a central galaxy or does not, while the num- 
ber of satellites is Poisson distributed about a mean A'sat- 
For each sample, we parameterize N{M) = Ncen + Naa.t with 
5 parameters 



iVcen(Af) = - erfc 



ln(Mcut/M) 



V2a 



and 



M - kMc, 



(Al) 



(A2) 



for M > wMcut and zero otherwise. Different functional 
forms have been proposed in the literature, but the current 
form is flexible enough for our purposes. 

Central galaxies always live at the minimum of the halo 
potential while satelli te galaxies are randomly placed assum- 
ing an NFW profile /Nava rro. Frenk fc Whit3 ll996'). If we 
instead use randomly chosen dark matter particles within 
halos, which preserves the anisotropy of the halos and any 
substructure, ^(r) is altered at the 10 per cent level on Mpc 
scales. The differences on large scales are very small. The 
concentration of the halo is taken from the N-body simu- 
lation, but multiplied by a free (mass independent) factor 
to allow galaxies to be more or less concentrated than the 
dark matter. This detail only affects the predictions on small 
scales. 

Given the 3D galaxy positions the correlation function 
is computed out to 10/i~^Mpc by direct pair counts, and 
then extrapolated assuming constant bias and a (dimension- 
less) mass power spectrum given by the Q-model 



A^(fc) = AL(fe) 



1 + 



(A3) 



1 + Ak 

with = 10 [2/(1 + z)]" '^^ and A = 1.7 ^"^Mpc. This form 
provides a reasonable fit to the DM power spectrum in the 
simulation over the redshift range of interest. The correla- 
tion function is then integrated, making the Limber approx- 
imation, to find u!{8s). 

Comparison of different paramerizations for Accn and 
Asat, different methods for making mock catalogs, different 
techniques for computing iL>{9s), different ranges and sub- 
sets of the data and the different simulations of the same 
cosmology, indicate that our results for the large-scale bias 
have systematic uncertainties at the several percent level. 
In the 500 /i"'^ Mpc box different realizations of the same 
HOD cause 2 per cent changes in the inferred large-scale 
bias, because of fluctuations in the galaxy-weighted mean 
halo mass from Poisson fluctuations in Asat. There is neg- 
ligible difference in the l/i~^Gpc box. By running a se- 
quence of boxes of increasing resolution, but with the same 
large-scale phases, we find that the cumulative mass func- 
tion changes by 5-10 per cent in the mass range of interest 



Slice 




ail 


(T 


K, 




X 


1 


13.15 


13.71 


1.16 


1.59 


0.81 


37.62 


2 


13.20 


14.00 


1.10 


0.37 


1.07 


26.75 


3 


13.06 


14.03 


0.13 


1.19 


1.38 


32.80 


4 


12.95 


14.05 


0.17 


0.64 


1.13 


33.25 


5 


12.98 


13.86 


0.85 


1.53 


0.96 


41.45 


6 


13.21 


13.90 


1.26 


1.18 


1.33 


62.41 



Table Al. The HOD parameters from our best fitting model 
and used to make the mock LRG catalogs. The meaning of the 
parameters is given in the text. The a-re for 25 data points. 



due to structural changes in the simulated halos with in- 
creasing force resolution. We might expect a similar change 
if we included baryonic cooling and star formation in our 
simulations. This affects the inferred number density for a 
given set of HOD parameters. However the positions, and 
hence clustering properites, of the halos are largely unaf- 
fected by increasing force resolution. In contrast, the bias of 
a given halo population is sensitive to finite box size effects 
on scales a few percent of the box size, and the sensitivity is 
larger the more biased the halo population under considera- 
tion. We choose to measure ^(r) in the 500 h~^Mpc box only 
out to 10/i~^Mpc because we find systematic differences in 
halo clustering between the large and small boxes for the 
rarer halos. This can be traced to the particular modes cho- 
sen in the initial conditions. If we restrict to 10 h~^Mpc the 
large-scale bias agrees between the two simulations to 2 per 
cent, less than the random error from the fits. The aver- 
age halo parameters also agree to within the chain-inferred 
dispersion. Unfortunately the halo bias is still slightly scale 
dependent at 10/i~^Mpc, as determined from our l/i~^Cpc 
simulation, so our results extrapolated assuming constant 
bias tend to over-estimate 6 by 5-10 per cent. We correct for 
this over-estimate for the values quoted in Tabled 

An investigation of all of these effects leads us to assign 
a 5 per cent systematic error bar to the large-scale LRG bias 
estimates we derive. This uncertainty, while comparable to 
or larger than the statistical error for the LRG sample, is 
irrelevant for our main conclusions. However, future work on 
modeling LRGs for galaxy formation and evolution will be 
limited by theoretical uncertainties, and not observational 
errors. For making mock catalogs, including QSOs, we use 
the 500/i~^Mpc box. This allows us to probe further down 
the mass function. The theoretical inaccuracies, of concern 
for the LRGs, are much smaller than the observational errors 
on the QSO-LRG cross-correlations. As shown in Figure |4j 
the mock catalogs produced in the 500 h~^Mpc box provide 
a very reasonable description of the LRG clustering on large 
and small scales. 



APPENDIX B: COUNTS OF NEIGHBOURS 

Previous authors (jSerber et al.l l2006l : 

IStrand. Brunner fc Mverd l2007l l estimated the clustering 
of quasars by measuring the overdensities of photometric 
galaxies around quasars in cylindrical apertures. While 
these measurements can be related to the correlation 
functions presented in this paper, there are a number 
of subtleties in their interpretation that have resulted in 
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confusion in the literature. This appendix attempts to 
clarify these measurements, as well as compare them with 

our results. 

Following ISerber et al.l (|2006l ). we define Nq^g,r{R) as 
the average number of photometric galaxies within cylindri- 
cal apertures of transverse physical (as opposed to angular) 
radii R, centered on quasars (q), spectroscopic galaxies (g) 
and random points (r). If we assume the photometric sam- 
ple has a normalized redshift distribution /(x), we use the 
formalism of 331 to write 



(Bl) 

(B2) 
(B3) 



Nr{R) = nrfip / 2'KR'dR! 
Jo 

and 

Nq,giR) = nq^gUp I 2tx R!dR' x 
Jo 

{R')\ , 

where nr,q.g is the number of random points, quasars and 
spectroscopic galaxies, rip is the areal density of the photo- 
metric sample, and w is the projected cross-correlation be- 
tween the different samples. As in ^ (fix)) implies the red- 
shift distribution of the photometric sample averaged over 
the redshift distribution of the spectroscopic sample. Nor- 
malizing the random galaxies to the number of spectroscopic 
targets, we then obtain 



Ng 

Nr 



= 1 + 



(fix)). 



2nR'dR' Wgp{R') 



(B4) 



for the quasars with a similar expression for galaxies. We 
theref ore s ee that the overdensities measured in lSerber et al.l 
l|2006l ) and lStrand. Brunner fc Mverj l|2007f ) can be related 
to the area-averaged projected correlation function weighted 
by the redshift distribution of the photometric galaxies at 
the redshifts of the spectroscopic targets. This implies that 
the overdensities of the quasars and galaxies cannot be di- 
rectly compared since they are scaled by different weights. 
Furthermore, the overall amplitude of the overdensities can- 
not be interpreted without knowledge of the redshift distri- 
bution of the photometric galaxies (which is non-trivial for 
high redshift quasars since one is starting to probe the high 
redshift tail of the photometric galaxies). 

Eq. IB4I has the property that it asymptotes to 1 on 
large scales irrespective of galaxy type. We estimate the 
scale at which Nq^g/Nr ~ 1 as follows - the redshift dis- 
tribution of the photometric galaxies spans approximately 
1 h~^Gpc and the integral of / must equal unity, suggesting 
(fix)) ~ 10~^ h Mpc~^ . Assuming that the cross-correlation 
between quasars (spectroscopic galaxies) and photometric 
galaxies is similar to the quasar-LRG cross-correlation, this 
implies that the second term of Eq. IB4I is 0(1) on scales 
R ~ 0.1/i~^Mpc. Note that this implies that the overdensi- 
ties are 0(1) on all scales, especially on scales larger than a 
Mpc. We emphasize that this is simply due to the division 
by the mean density and in no way implies that the quasars 
and spectroscopic galaxies have the same large scale bias or 
inhabit halos of similar masses. Note that this is a significant 
difference between spectroscopic and photometric samples - 
for spectroscopic samples, the width of the redshift distribu- 
tion is typically a few tens of Mpc (to integrate out redshift 
space distortions), and the second term in Eg. IB4l is much 




Figure Bl. The redshift distributions of spectroscopic L^, galax- 
ies [dotted, red], photome tric galaxies 19 < r < 21, [solid, black] 
jMandelbaum et al.l2007l ) and spectroscopic quasars with z < 0.4 
[dashed, blue]. Note that the normalizations of the redshift dis- 
tributions are arbitrary. The figure emphasizes the fact that 
galaxies and quasars probe very difi'erent regions of the photomet- 
ric redshift distribution, complicating the comparison between 
their clustering. 



larger than the first. In this case, the overdensities can be di- 
rectly interpreted as the angle- averaged correlation function. 
For photometric samples, the complications can be simply 
circumvented by subtracting 1 from the overdensities, if the 
redshift distribution of the photometric sa mple is known (see 
the discussion above). Howeve r , neit her iSerber et al.l (|2006l ) 
nor lStrand. Brunner &: IVIversI (|2007l ') do this when compar- 
ing with the galaxy samples, and their results must not be 
interpreted as quasars having the same clustering. An impor- 
tant corollary to this is that the upturn seen in overdensities 
cannot be interpreted (as has often been in the literature) 
as an excess in small-scale clustering, but is simply the sig- 
nature of a clustered sample of objects. Indeed, as we show 
below, the cross-correlation for the quasars is consistent with 
being a power-law down to small scales. 

The final complication in interpreting these results 
arises due to the width of redshift distribution of the photo- 
metric sample. Since the photometric sample covers a large 
redshift range, it cannot be treated as a homogeneous sam- 
ple and we cannot model its auto-correlation as we did with 
the LRGs. In order to interpret the cross-correlation of the 
quasars with the photometric sample, we must therefore 
compare it with the cross-correlation between a particular 
population of galaxies and the photometric sample. Unfortu- 
nately, for the SDSS, quasars and galaxies occupy different 
redshift ranges (see Fig. IBip and therefore probe different 
sub-populations of the photometric sample. Any comparison 
between quasars and galaxies must also take into account 
these population differences, complicating any analysis. 

A detailed modeling of all these effects goes well beyond 
the scope of this A ppendix. However, as an ill u strati on, we 
take the results of IStrand. Brunner fc IVIver^ (|2007f ). sub- 



© 0000 RAS, MNRAS 000, 000-000 



The real-space clustering of luminous red galaxies around z < 0.6 quasars in the Sloan Digital Sky Survey 



1000 


1 1 1 1 






' ' 1 


- 














(Mpc/ 












a 

1^ 

100 


', , , 1 






, , 1 





R (Mpc/h) 

Figure B2. The integrated, projected cross-correlation func- 
tion, ii)gp( < R), of photometric g alaxie s and quasars in 
SDSS from IStrand. Brunner fc Mveral ||2007| ') assuming f{xq) ~ 
10~'^hMpc~^ (see text). The quoted errors are the size of the 
plotted symbols. The line is derived from a power law correlation 
function with slope 1.9 and rg = 5 h~^Mpc. 

tract 1, and divide by {/(x)) ~ lO^'^/iMpc"^; the results 
are in Fig. IB2I The clustering is consistent with a three di- 
mensional power law of slope 1.9 and ro ~ 5/i~^Mpc. We 
remind the reader that the value of ro is completely degen- 
erate with our assumption for (fix)}- 
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